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Magnetic molecules adsorbed on a superconductor give rise to a local competition of Cooper 
pair and Kondo singlet formation inducing subgap bound states. For manganese-phthalocyanine 
molecules on a Pb(lll) substrate, featuring numerous different Kondo scales, such states are resolved 
by scanning tunneling spectroscopy. We show that numerical renormalization group calculations for 
an effective one channel Anderson impurity model explain the energy and weight of the bound states 
very well. The application of the model and its parameters are justified by scaling arguments. 
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Magnetic atoms and molecules adsorbed on supercon- 
ducting surfaces allow for the microscopic observation of 
the fascinating interplay of Cooper pair and Kondo sin- 
glet formation [THTOj. The ground state (GS) of such 
a combined system has been predicted to be either a 
Kondo screened singlet state (5 = 0), if the Kondo scale 
&b Tk is much larger than the superconducting pairing 
energy A sc , or an unscreened multiplet state (S > 0) 
for &b Tk <C A sc . Characteristic features include subgap 
bound states (BS) - often called Shiba states. For the sin- 
glet GS with screened impurity spin, the BS are S > 
excitations, where the Kondo singlet is broken, and for 
the multiplet GS the BS are singlet excitations includ- 
ing Kondo screening. When the energy scale for Kondo 
singlet formation becomes smaller and comparable with 
the Cooper pairing energy, fee Tk ~ A sc , the BS energy 
^ 0, and GS and BS cross. At this point a quan- 
tum phase transition from S = to S > occurs. The 
scattering and Cooper pair breaking effect is expected to 
behave like 1 — (i^/A sc ) 2 and thus most effective here 
leading to a strong suppression of the superconducting 
T c for larger impurity concentrations [4j [5j [7] . 

For a long time the accurate resolution of the sub- 
gap BS and their dependence on Tk and A sc have re- 
mained elusive. Recently, the BS have been analyzed in 
tunable mesoscopic superconductor-quantum dot-normal 
lead structures Scanning tunneling microscopy 

(STM) has been used to detect the local influence of 
single magnetic atoms on a superconducting substrate 
|12j . To study the interplay, an experiment with variable 
magnetic interaction strengths is desirable. Manganese- 
phthalocyanine (MnPc) molecules on a Pb(lll) substrate 
forms a Moire-like superstructure. Due to a large number 
of different adsorption sites, which lead to varying Kondo 
coupling strengths J, one can study different values of Tk 
and access the regime &b Tk ~ A sc [13J. Tunneling spec- 
troscopy on different Mn sites reveals varying and the 
occurrence of the predicted GS transition |13j . 



In spite of the complexity of MnPc the observations 
could be understood qualitatively by the physics of 
quantum impurity models [T3]. However, the question 
whether the experimental data agrees quantitatively with 
predictions from an effective one channel Anderson im- 
purity model (AIM), which has been subject of the the- 
oretical literature for decades, has remained open. To 
provide an answer is a crucial step to unite theory and 
experiment, and to validate the microscopic understand- 
ing of the physical processes of screening and pair break- 
ing. Relevant quantities to compare include the value 
^k,sc = ^B^k/Asc f° r the occurrence of the GS transi- 
tion and the positions E^ and weights of the BS. In 
this work we provide a detailed comparison of the ex- 
perimental observations with numerical renormalization 
group (NRG) [T4HI8] calculations, which are known to 
capture the Kondo effect accurately. Many other meth- 
ods have been used to describe theoretically impurities in 
superconductors on a microscopic scale, however, most of 
them contain mean field aspects like classical spins, such 
that the Kondo effect is not treated reliably pHSJ fl9l 
|2T] . For an appropriate parameter regime, we find good 
agreement with the experimental data, which supports 
the description in terms of an effective one channel model. 
Furthermore, we derive the effective model from the more 
general situation using scaling arguments. 

MnPc molecules (see Fig.[TJa)) have been deposited on 
an atomically clean Pb(lll) substrate at room tempera- 
ture under ultra-high vacuum conditions. STM at 4.5 K 
resolves highly ordered islands (Fig. [ljb)). Tunneling 
spectroscopy has been used to resolve the superconduct- 
ing gap structure and its subgap states, as well as Kondo 
resonances on different molecules. The superconducting 
state of the Pb(lll) substrate and Pb tip shows as pro- 
nounced differential conductance peaks at E = ±2A SC 
with A sc = l.lmeV. 

On the Mn sites, we identify a broad peak around the 
Fermi energy (see Fig. f3^b)) , which can be fitted by two 
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FIG. 1: (Color online) a) Molecular structure of Manganese- 
Phthalocyanine (MnPc), b) Constant-current STM image 
of highly ordered island of MnPc on Pb(lll) (I = 23pA, 
U=100mV), c) Tunneling spectra taken on different 
molecules on top of the Mn center (feedback opened at 
I = 450 pA, U = 8.6 mV) . Also shown are the fits used to derive 
the deconvoluted spectra below in (d) . e) Examples of NRG 
spectra with broadened peaks for similar Eb. 

Fano lineshapes, representing two different Kondo screen- 
ing processes with Tk,2 ^> ^k,i [13J. For the isolated 
MnPc complex, density functional theory (DFT) calcu- 
lations have found that the high spin configuration of Mn 
(S = 5/2) is reduced to S = 3/2 and the unpaired elec- 
trons occupy the b2 g (d xy ), the e g (d 7r )^ and the ai g (d z 2) 
orbital (x-axis along an arm of MnPc) and are aligned 
due to Hund's rule coupling [22J. When adsorbed on 
the Pb surface (in z direction), the d z i orbital hybridizes 
strongly with the Pb states and is therefore quenched 
|23| . Hence, two screening channels and a spin S = 1 
is left. Since &B^k,2 ^ ^sc? Kondo screening domi- 
nates for this channel |7J. Therefore, we only correlate 
Tk,i = Tk ~ A sc with the appearance of the BS. 

Three differential conductance spectra taken on differ- 
ent MnPc molecules are shown in Fig. [TJc). Two larger 
peaks as well as the two smaller peaks are located at 
symmetric bias voltages, while their intensity shows a 
strong asymmetry. Similar to the zero bias peak in the 
plain superconductor-superconductor junction, the small 
peaks are a result of thermal excitations across the gap. 
In order to remove the effect of the superconducting tip 
on the tunneling spectra, we developed a deconvolution 
method [13]. The result is representative for the quasi- 



particle density of states (DOS) of the MnPc molecule 
on the superconducting Pb surface (Fig. [TJd)) . These 
plots serve to deduce E b and their intensity. We observe 
a gradual increase in the asymmetry of the weights when 
the states shift closer to the Fermi level £p. This is shown 
in Fig.|2]as ratio w b ^ g /wb jS , where w b ^ g is the larger weight 
and w b ^ s the smaller one weight. 
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FIG. 2: (Color online) Ratio of bound state weights Wb, g /wb, s 
vs bound state energy Eb/A sc . Comparison of experimental 
results with theory. In all the range of bound state energy we 
used the model parameters £d/A sc ~ —73, and U/A sc ~ 91 
and the overlap F was varied; we did not consider any further 
variation of £d,U which could lead to a better agreement in 
some cases. 

Fig. [3] shows the dependence of E b on Rk,sc At large 
Rk,sc, the Kondo screening is efficient and the many- 
body ground state is a singlet. Here tunneling can oc- 
cur via the doublet state which involves breaking of the 
Kondo singlet and rearrangement of Cooper pairs. With 
decreasing Rk,sc this requires less energy and we find 
E b 0, and the BS and GS cross for R^ sc ~ 1.2. This 
point can be regarded as the critical point of the quan- 
tum phase transition and is a universal feature. Further 
decrease of Rk,sc leads to the unscreened doublet GS. 

For the theoretical description we use an AIM Hamil- 
tonian of the form 



H = H sc + Hd + H Y 
The superconducting medium reads, 



k,m,a 



k , rn 



(i) 

h.c.) 

where c\, m a creates a band electron with momentum 
fc, spin a and band index m, where m = 1, . . . N, and 
there are N available channels. Sk,m is the corresponding 
electronic dispersion and A sc the gap parameter chosen 
real. The band electrons hybridize with the impurity 
states via 
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FIG. 3: (Color online) a) Tunneling spectrum showing the 
Shiba states (feedback opened at I = 450pA, U = 8.6mV). b) 
Tunneling spectrum on the same molecule as in (a) in a larger 
energy scale evidencing a Kondo resonance (feedback opened 
at I = 470pA, U=130mV). c) Bound state energies E b /A sc 
vs Kondo temperature Tk/A sc . Comparison of experimental 
results with theory for the same choice of model parameters 
as in Fig. [2] 



where c\ G creates a d-level impurity electron with spin 
<j and index m. We will assume different matrix elements 
Vm due to different overlapping integrals. We define the 
width T m = irVmPm,c as the energy scale for the hy- 
bridization, p m?c is the DOS of the conduction band at 
£f- We will now show that the experimental results can 
be well described by a single channel model (N = 1). In 
this case the "d-orbital" term Hd simply reads, 



H d = ^2e d n a + Un^n^ 



(4) 



with the level position Ed relative to Sp and the on-site 
interaction with strength U, where n a = c^c^a- We 
set £y = in the following. The superconducting gap 
A sc = l.lmeV is fixed and sets the energy scale in all 
calculations. Even though the charge degrees of freedom 
of the impurity do not play a role at low energy, we stick 
to the description in terms of an AIM. The asymmetry 
of the weights can then be captured by a deviation 
from particle- hole symmetry, £ = Sd/U + 1/2 ^ 0. Alter- 
natively one could use a Kondo model with a potential 
scattering term. 

The choice of model parameters is motivated by the 
expectation that the main difference for the MnPc 
molecules in the different positions is the magnitude of 



the hybridization T, which in turn leads to different J 
and Tk- The energy Sd and the asymmetry £ are, on 
the contrary, expected to change little with the site [24] . 
Calculations for a set of parameters £^/A sc w —73, and 
U/A sc w 91, and varying r/A sc from 16 — 46 lead to an 
appropriate variation of and in Fig.[2]and give suit- 
able Kondo scales in Fig. [3j Good agreement of theory 
and experiment is obtained both for universal features 
such as i? K sc and non-universal ones like w^^/wb^. No- 
tice that Wb jg /wb jS increases whilst £ is constant. The 
smallest value for T corresponds to the lowest value of 
Tk, at the 5=1/2 doublet GS. The energy asymmetry 
of the levels is with £ ^ 0.3 quite large. Other methods 
[T91I2T] can give asymmetric weights even in more sym- 
metric situations, however, they are based on classical 
spins, such that the Kondo effect is not captured. 

Notice that the magnitudes of Sd, U and T do not 
correspond to usual atomic values ~ O(eV), but rather 
to values of the order lOOmeV. This is related to the fact 
that the AIM under consideration is an effective model 
valid for low energies (see discussion later). Calculations 
with Ed, U ~ 0(eV) produce a less pronounced increase 
of Wb^g/wb^ and can not explain the data very well. 

Of particular interest is the occurrence of the quantum 
phase transition. At first sight the experimental result for 
this transition, R^ sc — 1.2, in Fig. |3|seems at odds with 
earlier theoretical predictions. For ooth the one channel 
non- degenerate AIM and the Kondo model, NRG studies 
[T5UT8] have estimated that it occurs when sc ~ 0.3. 
We find that the origin of this discrepancy is the use of 
different definitions of Tk, which can vary by a prefactor. 
In the theoretical works [15-18J the definition [25J 



T K = 0.29[/7r] 1 / 2 e 



rre d (e d + U) 
2FU 



(5) 



was used. In this work, we employ the definition of Tk, 
widely used in experiments (26j[27], based on the width of 
the Kondo resonance Ak (half width at half maximum) 
in the limit T — >> 0. We adopt the same definition in 
our NRG calculations (Tk = Ak) and find that Ak and 
the definition in Eq. ([5| can differ by a factor 4 (see also 
Ref. |28j). Taking this into account our result for the 
quantum phase transition in Fig. [3] is in excellent agree- 
ment with the theoretical prediction for a one channel 
model. If more channels contribute for a given value of 
Tk, the transition occurs at smaller values of i? K sc [29] 
at variance with experiment. 

We now discuss the emergence of the low energy ef- 
fective one-channel model and the scaling of the model 
parameters. Using the DFT calculations for MnPc [52] 
as an input, we start with a model of the form of equation 
([I]) with N=2 and Ti<r2. Due to the spatial orientation 
the overlap is larger for the d n orbital. The impurity 
term can be written in terms of the level positions £d,m 
relative to £p, intra orbital U m , inter orbital Coulomb 
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FIG. 4: (Color online) Left: Schematic plot of the bare 
band and impurity energy scales and depiction of two scaling 
regimes (I, II). Right: Schematic plot of the scaling of the ef- 
fective level £d,i and Sd,i + Ui in the two regimes A £ (Di, D u ) 
and A £ (Tk,2, A)- Scaling stops at A ~ Jk,2, and the ef- 
fective one band model is used to compute the low energy 
properties. 

interaction U\2 and a Hund's rule interaction Jh, 

cr,ra m 

+ ^ Ul2nd J l J a^d J 2,a f ~ JftSd,! ' Sd,2, 



where S d,m 



y 



cr 2 ^d,m,cr 1 a cri,cr2^ d i' m i a 2 



Cd,i 



with the Pauli 



matrix a". The bare parameters of this model are not 
known, either from theoretical calculations or from exper- 
imental observations. However, generally we must have 
£d,m + U12 < and Sd y m + ^12 + U m > such that single 
occupation is favored. 

We now derive poor man's scaling equations [30] for 
the occupied level Sd,m and the interaction U m [25J. We 
focus on the quantities in channel 1. We only consider 
processes ~ V± explicitly, where the level occupations of 
6d,i changes. There are also contributions ~ V 2 2 , whose 
inclusion leads to quantitative changes in the equations 
but does not alter the conclusions. Band structure cal- 
culations [31] for Pb suggest a situation, where the band 
(—A, Ax) is asymmetric with respect to £p> A < Ax- 
Then, we have two scaling regimes (I, II, see Fig. ffl, and 
the scaling equations (I) for A G (A? At) read [9j|25], 



dsd,i 
dA 
dUi 



7T A-£ d> l' 

2r x c/i 



7T (A-e dj i)(A-e dj i)' 



(7) 
(8) 



where ^,1 = + C/i - Different from the usual ap- 
proaches we also use a scaling equation for T\, which is 
derived from dT^i/dA = 0, assuming Tk,i as a scaling 
invariant. Such a scaling procedure can be continued as 
long as the levels e^i, £d,i + U\ lie within (—A, A) and 
do not interfere with the Kondo scale. For \sd,i\ < A, 

de d ,i 



Ui\ < A, 



find 



> 0. This implies that the 



level €d,i is decreasing when scaling to lower energy. We 



also have ^ 



> and 



dA 



> 0, such that (e^i+L^i) 



decreases under the scaling. Hence, the model becomes 
more asymmetric in this regime. 

Once A ~ Di scaling changes to the equations (II), 



d£d,i 
dA 

dUi 
~dA 



Ti 2e d> i(e d ,i+A) + ^i(3e d ,i-A) 



7T 



(A 2 
4Ti ^d,i(A 2 - 



■4i)( A + £d,i) 
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For starting values, £d,i(A) < 0, e d ,i(Di) + /7i(A) > 0, 
and |e dj i| < A, kd,i + < A, one has < 0. 

In this case 6d,i is increasing under the scaling. For the 
same effective starting values we have > 0. Hence, 
C/i and (e dj i + C/i) decrease further under the scaling. 
This scaling can be continued until we reach A ~ Jh- 
We have Tk,2 ~ 50meF, A ~ 3eV, thus for usual esti- 
mates for Jh we expect Tk,2 < Jh < A- At this scale 
the spins lock into the high spin 5 = 1 configuration. As 
shown in Ref. [32], this leads to a reduction of the mag- 
netic coupling J — > J/ (25) and the Kondo scale. In our 
approach it can be included in Ti, such that the scaling 
is slowed down [32J. Scaling can then be continued until 
A ~ J2,k- Then one part of the impurity spin becomes 
Kondo screened and decouples. At this scale we have an 
effective single band model with 5 = 1/2. According to 
the described scaling procedure, the parameters are sub- 
stantially reduced from the bare values and generically 
asymmetric. This justifies our choice of parameters in 
the NRG calculations. 

In conclusion, we provide a unified experimental and 
theoretical perspective of the microscopic interplay of 
Kondo screening and superconducting pairing for MnPc 
on lead, as manifested in the subgap bound states. We 
show that in spite of the complexity of MnPc an effective 
description based on the one channel AIM captures both 
universal aspects like sc and non- universal ones like 
the asymmetry of the weights very well. In the future 
it would be interesting to envisage situations, where the 
impurity spins can interact, such that different kinds of 
quantum phase transitions can occur. 
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